function d = deviationlog10( s1, s2, stconc1, stconc2, time )       
    pts = 7:7:time;
    d1 = deval(s1,pts);
    d2 = deval(s2,pts);
    
    stRNATotal1 = RNA2RNATotal(stconc1(1), stconc1(10), stconc1(9));
    stRNATotal2 = RNA2RNATotal(stconc2(1), stconc2(10), stconc2(9));
    
    m = 0;
    
    for i = 1:length(pts)
        d1cur = d1(:,i);
        d2cur = d2(:,i);
        d1total = RNA2RNATotal(d1cur(1), d1cur(10), d1cur(9));
        d2total = RNA2RNATotal(d2cur(1), d2cur(10), d2cur(9));
        m = max( abs( log10(d1total / stRNATotal1)- log10(d2total / stRNATotal2) ), m );
    end
    
    d = m;
end
